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Summary. — I review several central open questions concerning GRB Jets. 

• I discuss a new estimate of the beaming correction for the rate of GRBs 
~ 75 ± 20. 

• I discuss the universal structured jet (USJ) model and conclude that while 
jets might be structured they are less likely to be universal. 

• I discuss recent observations of a sideways expansion of a GRB afterglow and 
compare these with current hydrodynamics simulations of jet evolution. 

• I discuss the implications of resent outliers to the energy-angle relation. 

PACS 98.70.Rz ,95.85.Pw,95.30.Lz Gamma ray bursts, gamma- 

rays, Hydrodynamics . 



1. Introduction 

The realization that the relativistic outflow in GRBs is in the form of jets has direct 
implication to their energy budget and their rates. This has, in turn, further implications 
on the nature of the inner engines. For example prior to the understanding that GRBs 
are beamed, events such as GRB 990123 with an isotropic equivalent energy of more than 
10 54 erg were hard to explain. With beaming the energy output of this event is a "mere" 
10 51 erg and it is compatible with a simple stellar mass progenitor. 

Evidence of jetted GRBs arises from observations H] of the predicted achromatic 
breaks in the afterglow light curves 0] . Further support is given by the comparison of 
long term radio calorimetry with the energy of the prompt emission 0. The time of the 
jet break provides an estimate of the jet angle 0]: 

(1) 9 = QA&r & d{nlE k ^ 52 ) l 'hl% iys = O^radK^^) 1 / 6 ^, 

where th,days is the break time in days and Ek,iso$i is the "isotropic equivalent" kinetic 
energy in units of 10 52 ergs, while g 52 is the real kinetic energy in the jet i.e: Ek e 52 — 
(0 2 /2)E k 
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Frail et al. |Jj and Panaitescu and Kumar |7] have estimated the opening angles 9 for 
several GRBs with known redshifts. They find that the total gamma-ray energy release, 
when corrected using a beaming factor, /j 

(2) Ej — fb-E-y.i so — —^-Eji SO: 

is clustered. This energy-angle relation is commonly called the Frail relation. A precursor 
of this discovery was found already in 1999 by Sari, Piran & Halprn j3] who found that 
the two brightest bursts known at that time were beamed. Bloom, Frail & Kulkarni 
8 confirmed this clustering around ~ 1.3 x 10 51 erg on a larger sample. This result 
is remarkable as it involves two seemingly unrelated quantities, 9 that is determined 
from the break in the afterglow light curve and E 1 ^ so which is a property of the prompt 
emission. The fact that the product of these two unrelated quantities is a constant is, in 
my mind, an indication that our overall model is correct. 

There are two leading models for the jet structure and for the interpretation of the 
jet break and the beaming angle. According to the original uniform jet model (UJ) the 
energy per solid angle is roughly constant within some finite opening angle, 9, and it 
sharply drops outside 9. Within the UJ model, the observed break corresponds to the 
jet opening angle, 9 (see Eq. 0. The Frail relation implies, here, that the total energy 
released in GRBs is constant and that the differences in the isotropic equivalent energies 
arise from variations in the opening angles. 

According to the alternative universal structured jet (US J) model all GRB 

jets are intrinsically identical. The energy per solid angle varies as a function of the angle 
from the jet axis. The jet break corresponds to the viewing angle of the observer and 
the Frail relation imposes a specific distribution of energy per unit angle in the jet: 



(3) £{9) = 



E /{it9 2 ), for9>9 c ; 
£ O /(7T0 C 2 ), for9<9 c , 



where Eq is a constant and 9 C is the core angle of the jet (TJ]|. While the UJ model 
contains a free function, the luminosity function or the corresponding opening angle 
distribution, the USJ model is complectly determined by the Frail relation. Its apparent 
luminosity function has the form &(L) oc L~ 2 We can use this lack of freedom 

within the USJ model to test it. 

The question which model is the correct one is still an open one. There are many 
other open questions concerning the structure of the jets and their evolution. In this talk 
I review some of these open questions. I refer the reader to a recent review |13j for a 
more extended overview on GRBs in general and on GRB jets in particular. 



2. The Beaming and the Rate of GRBs 

The overall GRB rate depends clearly on the amount that GRBs are beamed. Within 
the UJ model this has been measured traditionally in terms of the beaming correction 
factor, fr , which is defined as the ratio of total number of bursts to the observed ones. 
To estimate the overall GRB rate we need the average beaming correction such 
that ritrue = {fb 1 )nobs- The average is performed over the observed distribution. Taking 
into account the fact that for every observed burst there are / b _1 that are not observed, 
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Frail et al., estimated the average beaming correction^) 

(4) (A _1 )™ = ^E^- 520 ±85 > 

i 2 

where the sum is over the observed distribution. 

However, this calculation overestimates the actual beaming correction. In the intrinsic 
luminosity distribution there are many low luminosity bursts that have large opening 
angles. These bursts dominate the rate estimate. However, they are rather weak and can 
be observed only to small distances. Hence they are under- represented in the observed 
distribution |14j . This effect can be taken into account in the following way. For a given 
burst with a luminosity L we define the volume from which such a burst can be detected: 

f z(L) 

(5) V L = / dz(dV/dz)R GRB (z)/[R GRB (0)(l + z)} 1 

Jo 

where Rgrb(z) is the comoving rate of GRBs and z m (L) is the maximal redshift from 
which a burst with a luminosity L can be detected. Similarly Vx — Vl=oo is the whole 
effective volume of the observable universe. The intrinsic beaming correction can be 
written as: 

This estimate is, of course, somewhat model dependent as it requires an assumption on 
Rgrb ■ Guetta, Piran & Waxman ^2] find that for several models in which GRBs follow 
the SFR (fr )int — 75 ± 25, about a factor of 8 smaller than the previous estimate of 
520 ± 85 that did not take this effect into account. 

Following ^1] we can also estimate the beaming correction for the rate of GRBs within 
the USJ model in the following way. The total flux of GRBs per year (or per any other 
unit of time) is an observed quantity obtained by summing over the observed distribution. 
A comparison of this total flux with the total energy emitted by a single burst can tell 
us directly the total number of bursts. Integrating over the energy distribution (Eq. |3J) 
we obtain the total energy that a burst with a USJ emits: 

f&c /'Umax 

(7) E vs , = 2[ / (E o /e 2 c ) 2 9d0 + / E Q 6- l d6] = E [l + 21og(0 max /0 c )], 

Jo Je a 

where 9 ma x is the maximal angle to which the jet extends. This immediately implies 
that the ratio of UJs (emitting each Eq) to USJs (emitting each Eq[1 + 21og(0 max /0 c )]) 
required to explain the observed flux is 

(8) Nuj/Nusj = [1 + 2 log(6 max /6 c )}. 



( 1 ) Frail et al., 6 estimate the observed beaming factor distribution as p bs{fb) = {fo/ fo) a+1 
and {fb/fo) l3+1 for /;, /o respectively. They find that a is poorly constraint by the data 
while p = — 2.77lg'3 4 an< ^ W(/o) = — 2.91^o'o6 an d obtain (/ ! J~ 1 }foi by integrating over this 
distribution. 
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The upper and lower limits of this integral arc uncertain but the logarithmic dependancc 
implies that the factor cannot be smaller than 2 or much larger than 5. This implies 
that the number of USJs required to produce the observed flux is about factor of 4 below 
the corresponding number of UJs. Hence, the average beaming correction for USJs is 
- 20 ±10. 

3. Universal Structured Jets 

One of the intriguing open questions concerning GRB jets is their angular distribution. 
The two leading models are the UJ and USJ discussed earlier. The differences between 
USJ and UJ have crucial implications to the question of the nature of GRBs' inner 
engines and their progenitors. First, the universality of the USJ requires that more or 
less the same process operates with the same parameters within different GRBs. Second, 
a USJ carries roughly five times more energy than a UJ. This implies, for a USJ, an 
energy budget of ~ 10 52 erg. The rate of USJs is, correspondingly, smaller by a factor of 
five. It is therefore, important to ask whether there are observations that can distinguish 
between the two models. 




6 (rad) 

Fig. 1. - (a): The 2D distribution density, dn(z, ln9)/dzd\n 9, of the GRB rate as a function of 
2 and ln# in the USJ model. The white contour lines confine the minimal area that contains 
1 o of the total probability. The circles denote 16 bursts with known z and 9 0. (b): A limited 
redshift range, 0.8 < z < 1.7 (containing 10 out of the 16 data points) in which both redshift 
selections effects and the sensitivity to the unknown GRB rate are minimized. From |15| . 
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Perna, Sari & Frail J^l calculated the expected distribution of the observed opening 
angles within the USJ model assuming that GRBs follow the SFR. Following them we 
define the number of bursts in the interval (9, 9 + d9) and [z, z + dz) as: 

(9) —dddz = sl n9 RGRB{z) dV{z) T(9 z) 

( ' d9dz (1 + z) dz [ ' 1 

where, V{z) is the comoving volume element and T(9, z) depends on the distribution of 
GRB duration (see ^J^l f° r details). Perna et al., integrated over the redshift dis- 
tribution and found, with reasonable assumptions on T a remarkable agreement between 
the expected angular distribution dn(9)/d9 = J (dn/d9dz)dz and the observed angular 
distribution, lending a strong support to the USJ model. However, Nakar, Granot & 
Guetta compared of the two dimensional distribution dn/d9dz with the observed 
one. They found (see Fig. [3j that the two distribution disagree strongly. The observed 
points are very far from the location of the peak of the expected distributions. Some 
are even in a non-allowed region. This implies that the agreement between the observed 
angle distribution and the one predicted by the USJ model was just a coincidence and 
should not be taken as supporting this model. 

A similar discrepancy arose when Guetta et al., |14j compared the observed count 
(Cmax/Cmin) distribution for the BATSE long bursts sample with the one expected 
from the USJ model. Guetta et al., ^2 found that the USJ model predicts a significant 
excess of weak bursts as compared with the observed distribution. 

Selection effects and uncertainties mean that these discrepancies are insufficient to 
rule out the USJ model. Still these discrepancies certainly imply that the agreement 




Fig. 2. - The predicted differential C m ax/C m m distribution for the USJ model and the observed 
distribution taken from the BATSE catalog. The inconsistency at the low peak flux range is 
apparent. Predicted differential distributions for different parameterizations of the luminosity 
function within the UJ model are also shown. From |14) . 
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found for the angular distribution alone cannot be used, as hoped, to demonstrate the 
validity of this model. If it will turn out that these conclusions hold with additional data, 
we will have to reconsider the validity of the US J. Both discrepancies could be removed 
if we consider Structured Jets that are not universal. That is if we allow an angular 
dependence of the flow parameters (such as T and/or E) but we do not require that all 
jets are similar. Namely we could replace USJ with SJ. While such a solution is viable 
it clearly takes away some of the simplicity and the elegance of the USJ mode. 

4. Jet Evolution 

An important question that determines the observed light curve of the afterglow is 
the sideways expansion of the jet after T ~ where T is the Lorentz factor. The 
observations of the radio afterglow of GRB 030329 provided a unique opportunity to 
test this issue. Taylor et al., |20l I21| have measured the size of the radio afterglow of 
GRB 030329 between 20 and 300 days after the burst. Oren, Nakar & Piran [22] (see 
also [2S1 ) compared the observed sizes with several schematic models of spherical and 
jetted propagation. The remarkable results is that the apparent size of the afterglow is 
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Fig. 3. - upper panel: R± as a function of Td for different sideways expansion models in ISM. 



The energy to external density ratio E/n 
R 



0.6 ■ 10 erg cm and 6q = 0.06rad. Lower panel: 



as a function of Td for different opening angles in ISM, with a constant E/n 
3 . Tj is in days, j3 3 = ^thermal- From [55] 
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rather insensitive to the details of the model. This robustness is a good indication for the 
validity of the model, as it is difficult to force it to have different values and the values 
obtained fit the observations. On the other hand it is a drawback when one wishes to use 
the size to determine the parameters of the outflow. It is insensitive to these parameters. 

Oren et al., [221 find that the image of a spherically expanding fireball is largest with 
expansion as t 6 , while the size of a sideways expanding jet (at v w c) increases as t ' 5 . 
As can be seen in Fig. there is practically no difference between an expansion at the 
speed of light or at the sound speed, c-s/3. On the other hand the size of a non-expanding 
jet increases only as t Q/25 . A comparison of these models to the observations (see Fig. EJ) 
shows that the non-expanding jet model is inconsistent with the observations. While an 
addition of the faster expansion during the Newtonian phase |23| can somewhat alleviate 
the problem, it is not clear that this can be done with reasonable parameters |22| . 

Based on simple analytic model Sari Piran & Halpern 4 estimated that the relativis- 
tic jet will expand sideways almost at the speed of light. On the other hand Panaitescu 
& Meszaros \IQ estimated that the jet will expand only at the sound speed, c/ \/3- Both 
are consistent with the observations. More recently Kumar & Granot |TZ] integrated the 
hydrodynamic equations over the radial direction and obtained one dimensional simpli- 
fied hydrodynamics equations. Solving these equations they find no or little sideways 
expansion. Similar results were obtained in a two dimensional numerical integration of 
the full hydrodynamic equations |181 I19| . Remarkably this does not influence much the 
observed light curve (see |18p. 

We are left with a puzzle why do the numerical simulation indicate little or no expan- 
sion while the observations suggest a rather rapid sideways expansion. One may guess 
that the current computations are not refined enough to follow the jet evolution (see |2U 
for a detailed discussion of the problematics of these computations) . We probably have 
to wait for a higher resolution codes to resolve this problem. Another possibility is that 
the size of the radio afterglow does not trace well the size of the expanding jet. Detailed 
emission computations have to be carried out to determine this possibility. 

5. The Implication of Outliers to the Frail Relation and a Speculation 

We begun stressing the importance of the Frail (energy-angle) relation. However, 
Berger et al., (21] pointed out that in addition to the well known weak GRB 980425 there 
are three other outliers to this relation GRBs 980326, 980519 and 030329 . Intriguingly 
enough all outliers are weaker relative to the common value of ~ 10 51 erg. I would like 
to suggest a simple explanation to this phenomenon. It seems that the common value 
indicates indeed the available energy reservoir. Under optimal conditions a significant 
fraction of this energy is released as 7-rays. An intriguing question, incidentally, is how 
is this conversion so efficient? These efficient cases produce the brightest and easiest to 
detect bursts. The less efficient case are less powerful (in 7-rays) and hence are easily 
missed. GRB 980425 would not have been discovered if it was not so close. Another 
outlier, GRB 030329, is also nearby at z=0.168. The redshift of the other two outliers is 
unknown. Thus, I conclude with a speculation that as time progresses and more energies, 
redshifts and jet angles will become available it will turn out that the Frail relationship 
is satisfied as an inequality with ~ 10 51 erg being the upper limit to the 7-ray energy. 

* * * 

This work was supported by EU-RTN " GRBs Enigma and a Tool" and by US-Israel 
BSF. I thank J. Granot, D. Guetta, E. Nakar and E. Waxman for helpful discussions. 
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